SPECTRUM OF LARGE RANDOM REVERSIBLE MARKOV CHAINS: 
HEAVY TAILED WEIGHTS ON THE COMPLETE GRAPH 

CHARLES BORDENAVE, PIETRO CAPUTO, AND DJALIL CHAFAI 

Abstract. We consider the random reversible Markov kernel K obtained by assigning i.i.d. 
non negative weights to the edges of the complete graph over n vertices, and normalizing by 
the corresponding row sum. The weights are assumed to be in the domain of attraction of 
an o-stable law, a g (0,2). When 1 :C a < 2, we show that for a suitable regularly varying 
sequence k„ of index 1 — the limiting spectral distribution iia of KnK coincides with the 
one of the random symmetric matrix of the un-normalized weights (Levy matrix with i.i.d. 
entries). In contrast, when < a < 1, we show that the empirical spectral distribution of 
, K converges without rescaling to a non trivial law fla supported on [—1,1], whose moments 

are the return probabilities of the random walk on the Poisson weighted infinite tree (PWIT) 
pH ' introduced by Aldous. The limiting spectral distributions are given by the expected value of the 

' random spectral measure at the root of suitable self-adjoint operators defined on the PWIT. This 

characterization is used together with recursive relations on the tree to derive some properties 
of and Jia- We also study the limiting behavior of the invariant probability measure of K. 
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1. Introduction 

Let Gn = {Vm En) denote the complete graph with vertex set Vn = {!,..., n}, and edge 
set En — {{i,j}, 1 ^ i,j ^ n}, including loops 1 < i ^ n. Assign a non negative 

random weight (or conductance) Uij — Uj^i to each edge € En, and assume that the weights 

U = {Ui^; {h j} G En} are i.i.d. with common law £ independent of n. This defines a random 
network, or weighted graph, denoted (Gn, U). Next, consider the random walk on (G„, U) defined 
by the transition probabilities 

If. . " 
K,,j := , with p, := Yl U^-o- (1-1) 

The Markov kernel K is reversible with respect to the measure p = X^iev Pi^^ that 

for all i,j €E Vn- Note that we have not assumed that £ has no atom at 0. If = for some i, 
then for that index i we set Kij — 6ij, 1 ^ j ^ n. However, as soon as C is not concentrated 
at then almost surely, for all n sufficiently large, pi > for all 1 ^ i ^ n, is irreducible and 
aperiodic, and p is its unique invariant measure, up to normalization; see e.g. 

For any square nxn matrix M with eigenvalues Ai(M), . . . , A„(M), the Empirical Spectral 
Distribution (ESD) is the discrete probability measure with at most n atoms defined by 



1 " 



All matrices M to be considered in this work have real spectrum, and the eigenvalues will be 
labeled in such a way that A„(M) ^ • • • ^ Ai(Af). 

Note that K defines a square nxn random Markov matrix whose entries are not independent 
due the normalizing sums pi. By reversibility, K is self-adjoint in L^{p) and its spectrum cr(if) is 
real. Moreover, a{K) C [— and 1 G o-{K). Since K is Markov, its ESD pK carries further 
probabilistic content. Namely, for any £ G IN, iipe(i) denotes the probability that the random walk 
on (G„,U) started at i returns to i after £ steps, then the i^^ moment of pK satisfies 

r\'pK{dx) = -tT{K') = - Vp£(^). (1.2) 
J-i ^ " ^ 

Convergence of the ESD. The asymptotic behavior of pK as n — >■ oo depends strongly on the 
tail of C at infinity. When C has finite mean x C{dx) = m we set to = 1. This is no loss of 
generality since K is invariant under the dilation t — tUij. If C has a finite second moment we 
write (7^ = J^{x — 1)^ L{dx) for the variance. 

The following result, from [11 , states that if < ct^ < oo then the bulk of the spectrum of -JnK 
behaves, when n — > oo, as if we had truly i.i.d. entries (Wigner matrix). Without loss of generality, 
we assume that the weights U come from the truncation of a unique infinite table (C/ij)i,j^i of 
i.i.d. random variables of law C. This gives a meaning to the almost sure (a.s.) convergence of 
p^K- The symbol ^ denotes weak convergence of measures with respect to continuous bounded 
functions. Note that \i{^/nK) — yjn — ^ oo. 

Theorem 1.1 (Wigner-like behavior). //£ has variance < < oo then a.s. 

1 " 

PV^K - n ^ ^v^A,(K) ^ , (1.3) 

k=l 

where W2ct *s the Wigner semi-circle law on [— 2(T, +2(t]. Moreover, if C has finite fourth moment, 
then X2{\/nK) and Xn{^/nK) converge a.s. to the edge of the limiting support [—2a,+2a]. 

This Wigner-like scenario can be dramatically altered if we allow £ to have a heavy tail at 
infinity. For any a € (0, oo), we say that C belongs to the class if C is supported in [0, oo) and 
has a regularly varying tail of index a, that is for all t > 0, 

G(<) := C{{t, oo)) L{t) t-'^ (1.4) 
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where L is a function with slow variation at oo, i.e. for any a; > 0, 

hm ^ ^ 1. 

t-!-00 L{t) 

Set a„ = inf{a > : nG{a) ^ 1}. Then nG(a„) — nL{a„)a^°' — > 1 as n — !■ oo, and 

nG{ant) — as n —> oo for all i > 0. (1.5) 
It is well known that a„ has regular variation at oo with index 1/a, i.e. 

for some function I with slow variation at oo, see for instance Resnick [531 Section 2.2.1]. As an 
example, if V is uniformly distributed on the interval [0, 1] then for every a G (0, oo), the law of 
supported in [l,oo), belongs to Hq. In this case. Lit) = 1 for i > 1, and a„ — m}^". 
To understand the limiting behavior of the spectrum of K in the heavy tailed case it is important 
to consider first the symmetric i.i.d. matrix corresponding to the un-normalized weights Ui^j. More 
generally, we introduce the random n x n symmetric matrix X defined by 

X = {Xij)i^ij^n (1.6) 

where (^i jOi^i^j^n are i.i.d. such that Uij :~ \Xij\ has law in Hq with a £ (0,2), and 

e= lim-&^e[0,l]. (1.7) 

t-)-oo r[\Xij\ > t) 

It is well known that, for a € (0,2), a random variable Y is in the domain of attraction of an 
a-stable law iff the law of |y| is in Hq, and the limit (II. 7p exists, cf. |17[ Theorem IX. 8. la]. It 
will be useful to view the entries Xij in (jl.6p as the marks across edge {i,j} € En of a random 
network (G„,X), just as the marks Uij defined the network (G„,U) introduced above. 

Remarkable works have been devoted recently to the asymptotic behavior of the ESD of matrices 
X defined by (|1.6p . sometimes called Levy matrices. The analysis of the Limiting Spectral Dis- 
tribution (LSD) for a E (0,2) is considerably harder than the finite second moment case (Wigner 
matrices), and the LSD is non explicit. Theorem 11.21 below has been investigated by the physicists 
Bouchaud and Cizeau [15] and rigorously proved by Ben Arous and Guionnet [7:, and Belinschi, 
Dembo, and Guionnet [5]; see also Zakharevich for related results. 

Theorem 1.2 (Symmetric i.i.d. matrix, a € (0,2)). For every a £ (0,2), there exists a symmetric 
probability distribution fia on II depending only on a such that (with the notations of U.ML6]) a.s. 



1 " 



Ma- 

71 ^ — ^ --y n— >-oo 

1=1 

In Section 13.21 we give a new independent proof of Theorem 11.21 The key idea of our proof 
is to exhibit a limiting self-adjoint operator T for the sequence of matrices a~^X, defined on 
a suitable Hilbert space, and then use known spectral convergence theorems of operators. The 
limiting operator will be defined as the "adjacency matrix" of an infinite rooted tree with random 
edge weights, the so called Poisson weighted infinite tree (PWIT) introduced by Aldous [T], see 
also [3|. In other words, the PWIT will be shown to be the local weak limit of the random network 
(G„,X) when the edge marks Xij are rescaled by a„. In this setting the LSD fia arises as the 
expected value of the (random) spectral measure of the operator T at the root of the tree. The 
PWIT and the limiting operator T are defined in Section [51 Our method of proof can be seen as 
a variant of the resolvent method, based on local convergence of operators. It is also well suited 
to investigate properties of the LSD fia, cf. Theorem II .61 below. 

Let us now come back to our random reversible Markov kernel K defined by p.l[) from weights 
with law C G Ha- We obtain different limiting behavior in the two regimes a G (0, 1) and a S (1, 2). 
The case a > 2 corresponds to a Wigner-type behavior (special case of Theorem I l.ip . We set 

— riCL^ . 
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Theorem 1.3 (Reversible Markov matrix, a € (1, 2)). Let fj.^ be the probability distribution which 
appears as the LSD in the symmetric i.i.d. case (Theorem \1.2\l . If L £ Hq, with a G (1, 2) then a.s. 



1 " 



Theorem 1.4 (Reversible Markov matrix, a € (0, 1)). For every a € (0, 1), there exists a sym- 
metric probability distribution Jla supported on [—1, 1] depending only on a such that a.s. 

1 " 

n '—^ ■ ' ri— >oo 
fc=l 

The proof of Theorem 11.31 and Theorem 11.41 is given in Sections 13.31 and 13.11 respectively. As 
in the proof of Theorem 11.21 the main idea is to exploit convergence of our matrices to suitable 
operators defined on the PWIT. To understand the scaling in Theorem ll.3l we recall that if a > 1, 
then by the strong law of large numbers, we have n~^pi — 1 a.s. for every row sum pi, and this is 
shown to remove, in the limit n — !■ oo, all dependencies in the matrix na~^K, so that we obtain 
the same behavior of the i.i.d. matrix of Theorem 11.21 On the other hand, when a G (0, 1), both 
the sum pi and the maximum of its elements are on scale a„. The proof of Theorem II .41 shows that 
the matrix K converges (without rescaling) to a random stochastic self-adjoint operator K defined 
on the PWIT. The operator K can be described as the transition matrix of the simple random 
walk on the PWIT and is naturally linked to Poisson-Dirichlet random variables. This is based 
on the observation that the order statistics of any given row of the matrix K converges weakly to 
the Poisson-Dirichlet law PD(a, 0), see Lemma [2.41 below for the details. In fact, the operator K 
provides an interesting generalization of the Poisson-Dirichlet law. 

Since is supported in [—1, 1], (|1.2I) and Theorem 11.41 implv that for all £ ^ 1, a.s. 

-y^P<?(*)= [ x'^p.Kidx) — !• I x'''Jla{dx) =: . (1.8) 



i=l 



R 



The LSD pa will be obtained as the expectation of the (random) spectral measure of K at the root 
of the PWIT. It will follow that 7^ (the moment of pa) is the expected value of the (random) 
probability that the random walk returns to the root in ^-steps. In particular, the symmetry of pa 
follows from the bipartite nature of the PWIT. 

It was proved by Ben Arous and Guionnet [71 Remark 1.5] that a G (0, 2) 1— >■ pa is continuous 
with respect to weak convergence of probability measures, and by Belinschi, Dembo, and Guionnet 
[51 Remark 1.2 and Lemma 5.2] that pa tends to the Wigner semi-circle law as a 2. We believe 
that Theorem 11.31 should remain valid for a — 2 with LSD given by the Wigner semi-circle law. 
Further properties of the measures pa and pa are discussed below. 

The case a = 1 is qualitatively similar to the case a € (1, 2) with the difference that the sequence 
Kn in Theorem 11.31 has to be replaced by k„ = na~^Wn where 

Wn = I xC{dx) . (1-9) 
Jo 

Indeed, here the mean of Ui,j may be infinite and the closest one gets to a law of large numbers is 
the statement that pi/nwn — >■ 1 in probability (see Section l3.4p . The sequence Wn (and therefore 
K,n) is known to be slowly varying at 00 for a = 1 (see e.g. Feller [171 VIII. 8]). The following mild 
condition will be assumed: There exists < £ < 1/2 such that 

liminf^^>0. (1.10) 

n— >oo Wn 

For example, if is uniform on [0, 1], then k„ = Wn — logn and lim„_i.oo wl„sj /w„ = e. In the 
next theorem pi stands for the LSD pa from Theorem 11.21 at a = 1. 

Theorem 1.5 (Reversible Markov matrix, a = 1). Suppose that C G with a = 1 and assume 
(jl.lOp . If Pk,„k is the BSD of KnK, with Kn = na^^Wn, then, as n ^ 00, a.s. Pk.„k Mi ■ 
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Properties of the LSD. In Section |4] we prove some properties of the LSDs fia and ^a- 

Theorem 1.6 (Properties of /Iq). Let be the symmetric LSD in Theorems M.'MTT^ 

(i) /Iq is absolutely continuous on R with hounded density. 

(ii) The density of /Iq at is equal to 



TT V Vr(i + ^ 

(iii) fj,a is heavy-tailed, and as t goes to +oo, 

/ia((i,+00)) - V". 

Statements (i)-(ii) answer some questions raised in [7J|3]. Statement (iii) is already contained 
in [5j Theorem 1.7], but we provide a new proof based on a Tauberian theorem for the Cauchy- 
Stieltjes transform that may be of independent interest. 

Theorem 1.7 (Properties of pta). Let Jla be the symmetric LSD in Theorem \1.4\ with moments 
7f as in (|1.8p . Then the following statements hold true. 

(i) For a G (0, 1), there exists 6 > such that 

72n ^ 5n~°' for all n ^ 1 . 

Moreover, we have liminfa/'i 72 > 0. 

(ii) For the topology of the weak convergence, the map a t-^ fla is continuous in (0, 1). 

(iii) For the topology of the weak convergence, 

hm Ha = + -do + -Oi. 

a\Q 4 2 4 

It is delicate to provide liable numerical simulations of the ESDs. Nevertheless, Figure [Uprovides 
histograms for various values of a and a large value of n, illustrating Theorems 1 1 . 3l - [T77l 

Invariant measure and edge-behavior. Finally, we turn to the analysis of the invariant probability 
distribution p for the random walk on (G,U). This is obtained by normalizing the vector of row 
sums p: 

P== (pi H hp„)"^(Pi>---;Pn)- 

Following [111 Lemma 2.2], if a > 2 then nmaxi^i^„ \pi — n^^\ — > as 71 — > cxi a.s. This uniform 
strong law of large numbers does not hold in the heavy-tailed case a E (0, 2): the large n behavior 
of p is then dictated by the largest weights in the system. 

Below we use the notation p — (pi, . . . , p„) for the ranked values of pi, . . . , /5„, so that pi > 
P2 ^ • • • and their sum is 1. The symbol denotes convergence in distribution. We refer to 
Subsection 12.41 for more details on weak convergence in the space of ranked sequences and for the 
definition of the Poisson-Dirichlet law PD{a, 0). 

Theorem 1.8 (Invariant probability measure). Suppose that C G Hq,. 

(i) If a £ {0,1), then 

P ^ ^ iVl,VuV2,V2,...) , (1.11) 

where Vi > V2 > ■ ■ ■ stands for a Poisson-Dirichlet PD(a,0) random vector. 

(ii) Ifae (1,2), then 

Kn{n+1)/2P lr{xi,Xi,X2,X2,...) , (1.12) 

where Xi > X2 > ■ ■ ■ denote the ranked points of the Poisson point process on (0,oo) with 
intensity measure ax^^^^dx. Moreover, the same convergence holds for a = 1 provided 
the sequence k„ is replaced by na^^Wn, with Wn as in il.9\) . 

Theorem 11.81 is proved in Section \S\ These results will be derived from the statistics of the 
ranked values of the weights Uij, i < j, on the scale ara(n+i)/2 (diagonal weights Uij are easily 
seen to give negligible contributions). The duplication in the sequences in (jl.l2p and (jl.lll) then 
comes from the fact that each of the largest weights belongs to two distinct rows and determines 
alone the limiting value of the associated row sum. 
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Theorem 11.81 is another indication that the random walk with transition matrix K shares the 
features of a trap model. Loosely speaking, instead of being trapped at a vertex, as in the usual 
mean field trap models (see [TH [SI [THl [IB]) here the walker is trapped at an edge. 

Large edge weights are responsible for the large eigenvalues of K . This phenomenon is well 
understood in the case of symmetric random matrices with i.i.d. entries, where it is known that, 
for a E (0,4), the edge of the spectrum gives rise to a Poisson statistics, see [261 H]. The behavior 
of the extremal eigenvalues of K when C has finite fourth moment has been studied in In 
particular, it is shown there that the spectral gap 1 — A2 is 1 — 0{n~^/^). In the present case of 
heavy-tailed weights, in contrast, by localization on the largest edge-weight it is possible to prove 
that, a.s. and up to corrections with slow variation at 00 

i-A.=|^("";rw "'^'''^ (1-13) 

Similarly, for a € (2,4) one has that A2 is bounded below by ri~("~^^/". Understanding the 
statistics of the extremal eigenvalues remains an interesting open problem. 

2. Convergence to the Poisson Weighted Infinite Tree 

The aim of this section is to prove that the matrices X and K appearing in Theorems 11.21 11.31 
and 11.41 when properly rescaled, converge "locally" to a limiting operator defined on the Poisson 
weighted infinite tree (PWIT). The concept of local convergence of operators is defined below. We 
first recall the standard construction of the PWIT. 

2.1. The PWIT. Given a Radon measure ly on R, PWIT(i^) is the random rooted tree defined 
as follows. The vertex set of the tree is identified with IN-'' UfegnlN'^ by indexing the root as 
]N° = 0, the offsprings of the root as IN and, more generally, the offsprings of some v € IN'^' as 
(vl), (v2), • • • g IN'^+^ (for short notation, we write (vl) in place of (v, 1)). In this way the set of 
V G IN" identifies the n^^ generation. 

We now assign marks to the edges of the tree according to a collection {Sv}vgin/ of independent 
realizations of the Poisson point process with intensity measure v on R. Namely, starting from the 
root 0, let S0 = {yi, j/2, • • ■ } be ordered in such a way that | ^ |y2| =^ • • • , and assign the mark 
Ui to the offspring of the root labeled i. Now, recursively, at each vertex v of generation k, assign 
the mark y^i to the offspring labeled vi, where Sv = {j/vi, 2/v2, • ■ • } satisfy \yvi\ ^ |yv2| 

Note that Sv has in average r^(Il) £ (0,oo] elements. As a convention, if z^(]R) < 00, one sets 
the remaining marks equal to 00. For example, ii v = X5i is proportional to a Dirac mass, then, 
neglecting infinite marks, PWIT(i/) is the tree obtained from a Yule Process (with all marks equal 
to 1). In the sequel we shall only consider cases where ly is not finite and each vertex has a.s. an 
infinite number of offsprings with finite and distinct marks. If ly is the Lebesgue measure on [0, cxi) 
we obtain the original PWIT in [I]. 

2.2. Local operator convergence. We give a general formulation and later specialize to our 
setting. Let V he a, countable set, and let L'^{V) denote the Hilbert space defined by the scalar 
product 

(0, tp) := X! ' " '^^ 

where (f>,ip £ and Su denotes the unit vector with support u. Let V denote the dense subset of 
L^{V) of vectors with finite support. 

Definition 2.1 (Local convergence). Suppose S„ is a sequence of bounded operators on LF'{V) and 
S is a closed linear operator on Lp'iV) with dense domain -D(S) D 2?. Suppose further that D is 
a core for S (i.e. the closure of S restricted to V equals S ). For any u,v £ V we say that (S„, u) 
converges locally to {S,v), and write 

(S„,m) ->■ (S.v) , 

if there exists a sequence of bijections Un '■ V ^ V such that crn{v) ~ u and, for all cj) ED, 
in L'^{V), as n ^ 00. 
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In other words, this is the standard strong convergence of operators up to a re-indexing of V 
which preserves a distinguished element. With a shght abuse of notation we have used the same 
symbol (t„ for the linear isometry (T„ : L'^{V) — > L'^{V) induced in the obvious way, i.e. such that 
(TnSv = ^a„(v) all V €V. The point for introducing Definition 12.11 lies in the following theorem 
on strong resolvent convergence. Recall that if S is a self- adjoint operator its spectrum is real and 
for all z E C+ := {z G C : 3z > 0}, the operator S — z/ is invertible with bounded inverse. The 
operator-valued function z i— >■ (S — zl)~^ is the resolvent of S. 

Theorem 2.2 (From local convergence to resolvents). //S„ and S are self-adjoint operators that 
satisfy the conditions of Definition \2.1\ and (S„, m) (S, v) for some u,v eV , then, for all z e 

{5u. (S„ - ziyHu) ^ {Sv, (S - zI)-H.,). (2.1) 

Proof of Theorem \2.SX It is a special case of [23l Theorem VIII. 25(a)]. Indeed, if we define S„ = 
a~^SnO'n, then S„(/< — >■ S0 for all (/) in a common core of the self-adjoint operators S„,S. This 
implies the strong resolvent convergence, i.e. (S„ — zT)^^%Ij — >■ (S — zI)^^%Ij for any z G C+, 
■0 G L?'{V). The conclusion follows by taking the scalar product 

(<5„ (S„ - zI)-H,) = {5u, (S„ - zI)-Hu)- 

□ 

We shall apply the above theorem in cases where the operators S„ and S are random operators 
on L^(y), which satisfy with probability one the conditions of Definition 12. II In this cases we say 
that (S„,u) —7- (S,f) in distribution if there exists a random bijection cr„ as in Definition 12.11 such 
that (T~^S„cr„(/) converges in distribution to S0, for all G 2? (where a random vector -!/;„ G L'^iV) 
converges in distribution to ip if 

lim E/(i/-„) = KfiiP) 

n— f oo 

for all bounded continuous functions / : L'^{V) — R). Under these assumptions then (12. ip becomes 
convergence in distribution of (bounded) complex random variables. In our setting the Hilbcrt 
space will be L^{V), with V = W-^, the vertex set of the PWIT, the operator S„ will be a rescaled 
version of the matrix X defined by (II. 6p or the matrix K defined by (|l.ip . The operator S will be 
the corresponding limiting operator defined below. 

2.3. Limiting operators. Let 9 be as in Theorem 1 1.2 1 and let £g be the positive Borel measure on 
the real line defined by d£g{x) = dl^x-^Qjdx + {l — d)l^^^Q^dx. Consider a realization of PWIT(^e). 
As before the mark from vertex v G f^'^ to vfc G W'^^^ is denoted by j/v/c- We note that almost 
surely 

^|yvfcr'/"<cx3, (2.2) 

k 

since a.s. lim?; \yvk\/k — 1 and 'J2k^^^^°' converges for a G (0,2). Recall that for V = fi^ , V 
is the dense set of L'^{V) of vectors with finite support. We may a.s. define a linear operator 
T :V ^ L'^iV) by letting, for v, w G M-^', 

{sign(?/w)|j/w|~^^" if w = vfc for some integer k 
sign(?/v)|2/v|^^^" if V = wfc for some integer k (2-3) 
otherwise. 

Note that if every edge e in the tree with mark ye is given the "weight" sign(?/e)|ye| ^^^" then we may 
look at the operator T as the "adjacency matrix" of the weighted tree. Clearly, T is symmetric, 
and therefore it has a closed extension with domain D{T) C L'^{¥(^) such that T> C D{T); see 
e.g. [221 VIII. 2]. We will prove in Proposition IA.2I below that T is essentially self-adjoint, i.e. 
the closure of T is self-adjoint. With a slight abuse of notation, we identify T with its closed 
extension. As stated below, T is the weak local limit of the sequence of n x n i.i.d. matrices a~^X, 
where X is defined by (jl.6p . To this end we view the matrix X as an operator in L'^{V) by setting 
{Si, X6j) — Xij, where i, j G M denote the labels of the offsprings of the root (the first generation), 
with the convention that Xij — when either i > n or j > n, and by setting {S-a, XS^) — when 
cither u or v does not belong to the first generation. 
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Similarly, taking now = 1, in the case of Markov matrices K defined by (jl.lL for a e [1, 2), 
T is the local limit operator of KnK. To work directly with symmetric operators we introduce the 
symmetric matrix 

5 . = (2.4) 

which is easily seen to have the same spectrum of K (see e.g. Lemma 2.1]). Again the matrix 
S can be embedded in the infinite tree as described above for X . 

In the case a g (0, 1) the Markov matrix K has a different limiting object that is defined as 
follows. Consider a realization of PWIT(^i), where ^l is the Lebesgue measure on [0,oo). We 
define an operator corresponding to the random walk on this tree with conductance equal to the 
mark to the power —1/a. More precisely, for v G IN-'', let 

p(v) = + 

feeiN 

with the convention that y^^^" = 0. Since a.s. limfe |yvfe|/fc = 1, p(v) is almost surely finite for 
a & (0, 1). We define the linear operator K on V, by letting, for v, w e IN-'', 

s if w = vfc for some integer fc 

p(v) tJ 

K(v, w) = ((5v, K^w) = \ i^Y-r if v = wfc for some integer fc (2-5) 



P(v) 

otherwise. 



Note that K is not symmetric, but it becomes symmetric in the weighted Hilbert space L^{V,p) 
defined by the scalar product 



-2/ 



Moreover, on L {V,p), K is a bounded self-adjoint operator since Schwarz' inequality implies 

U V 



so that the operator norm of K is less than or equal to 1. To work with self-adjoint operators in 
the unweighted Hilbert space L'^{V) we shall actually consider the operator S defined by 



S(v,w):../44k(v,w) = -^^. (2.6) 



This defines a bounded self-adjoint operator in L^{V). Indeed, the map — > ^/p(y)Sv induces a 
linear isometry D : L^(V, p) — L'^iV) such that 

(0,S^) = (D-V,KD-V)p, (2.7) 

for all (/), ?A S L'^{V). In this way, when a G (0, 1), S will be the limiting operator associated to the 
matrix S defined in (|2.4p . Note that no rescaling is needed here. The main result of this section is 
the following 

Theorem 2.3 (Limiting operators). As n goes to infinity, in distribution, 

(i) ifae (0, 2) and 9 E [0, 1], then {a-^X, 1) (T, 0), 

(ii) ifa€ (1, 2) and 9 ^ 1 then (nnS, 1) (T, 0), 

(iii) ifae (0, 1) then {S, 1) ^ (S, 0). 

From the remark after Theorem l2.2l we see that Theorem 12 .31 implies convergence in distribution 
of the resolvent at the root. As we shall see in Section [3l this in turn gives convergence of the 
expected values of the Cauchy-Stieltjes transform of the BSD of our matrices. The rest of this 
section is devoted to the proof of Theorem [ 
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2.4. Weak convergence of a single row. In this paragraph, we recall some facts about the 
order statistics of the first row of the matrix X and K, i.e. 

. . . ,Xi^„) and (C/i,i, . . . , [/i,„)/pi , 

where Ui,j = \Xi^j\ has law Hq. Let us denote by ^ V2 • • • 14, the order statistics of the 
variables Uij, 1 < j < n. Recall that pi = X]j=i ^j- Let us define Afe^„ = X]j=fc+i ^ fo^' k < n 
and A| „ J2^=k+i ^j- ^he set of sequences {vj} S [0, 00)*' with vi ^ V2 ^ ■ ■ ■ ^ such 

that hmj_i.oo Vj — 0, and let ^1 C -4 be the subset of sequences satisfying J2j — 1- We shall 
view 

y„=(— ) and Z„= (—,...,— 

\an an J \Pi Pi 

as elements of A and Ai, respectively, simply by adding zeros to the right of Vn/un and Vn/pi- 
Equipped with the standard product metric, A and Ai are complete separable metric spaces (^1 
is compact) and convergence in distribution for A, y^i-valued random variables is equivalent to 
finite dimensional convergence, cf. e.g. Bertoin [9|. 

Let £'i,i?2,... denote i.i.d. exponential variables with mean 1 and write 7^ — We 
define the random variable in A 

The law of Y is the law of the ordered points of a Poisson process on (0, 00) with intensity measure 
ax^^^^dx. For a £ (0, 1) we define the variable in Ai 



Z = 



7i 72 



n=l 7" 2^ra=l 7" 



For a G (0, 1) the sum 7^ ° is a.s. finite. The law of Z in Ai is called the Poisson-Dirichlet 
law PD(a,0), see Pitman and Yor [351 Proposition 10]. The next result is rather standard but we 
give a simple proof for convenience. 

Lemma 2.4 (Poisson-Dirichlet laws and Poisson point processes). 

(i) For all a > 0, Yn converges in distribution to Y. Moreover, for a € (0,2), (a^^Pj)j^i is 
a.s. uniformly square integrable, i.e. a.s. limfcSup„>;. „ ~ 0. 

(ii) If a G (0,1), Zn converges in distribution to Z. Moreover, {o.n^Vj)j^i is a.s. uniformly 
integrable, i.e. a.s. lim^ sup„>j. Afe^„ =0. 

(iii) If I C K is a finite set and ^ V2 ^ • • • denote the order statistics 0/ {f/i.j}jG{i....,n}\/ 
then (i) and (ii) hold with Y^ = (V^//a„, V//a„, . . . ) and Zf^ — {yl j p\, Vl j p\, . . . ). 

As an example, from (i), we retrieve the well-known fact that for any a > 0, the random 
variable max(t/i.i, . . . , [/i,„) converges weakly as n — > 00 to the law of 7]^ " . This law, known 
as a Frechet law, has density ax~'^~^ on (0, 00). 



Proof of Lemma \2.4\ As in LePage, Woodroofe and Zinn [50] we take advantage of the following 
well known representation for the order statistics of i.i.d. random variables. Let G be the function 
in (|1.4p and write 

G-\u) ^ M{y > : G{y) ^ u} , 
u E (0, 1). We have that {Vi, . . . ,Vn) equals in distribution the vector 

(G-' (7i/7n+i) , . . . , (7«/7«+i)) , (2.8) 
where "fj has been defined above. To prove (i) we start from the distributional identity 

^ G-M7l/7n+l) G-Hlnhn+l Y 
— I 5 • ■ • 1 

which follows from (|2.8p . It suffices to prove that for every k, almost surely the first k terms above 
converge to the first k terms in Y. Thanks to (|1.5p . almost surely, for every j: 

a-^G-i(7,/7n+i)^77" , (2.9) 
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and the convergence in distribution of y„ to Y follows. Moreover, from (jl.Sp . for any (5 > we can 
find no such that 

— 1 —1 —1 — — 

o-n yj = 0'n G (7j/7„+i) < (n7j/(l + (5)7„+i) ° , 

for n 5^ no, 7 G IN. Since n/^n+i 1, a.s. we see that the expression above is a.s. bounded by 

2(1 + (5)a7^ ° , for n sufficiently large, and the second part of (i) follows from a.s. summability of 

__2 

Similarly, if a G (0, 1), Afe_„ has the same law of 

^ ilj/ln+l) , 

j=k+l 

_ j_ 

and the second part of (ii) follows from a.s. summability of 7^ " . To prove the convergence of Z„ 
we use the distributional identity 



7 - 



G-M7l/7n+l) G-l(7n/7n+l) 



t:u (7./7n+i) ' ■ ■ ■ ' (7,/7n+i) ; ■ 

As a consequence of (12. 9p . we then have almost surely 

n oo ^ 

<^ il'jhn+l) ^ 51 ~ ' 

i=i i=i 
and (ii) follows. Finally, (iii) is an easy consequence of the exchangeability of the variable [Ui^i) : 

+ Vk) ^ P(3j G / : f/ij ^ Vk) SC |/|P(C/i,i ^ 14) = |/|^. 

□ 

The intensity measure ax~'^~^dx on (0, oo) is not locally finite at 0. It will be more convenient 
to work with Radon (i.e. locally finite) intensity measures. 

Lemma 2.5 (Poisson Point Processes with Radon intensity measures). Let f", > • ■ • sequences 
of i.i.d. random variables on R := RU {±00} such that 



nP(CrG-) ^ (2.10) 
where v is a Radon measure on R. Then, for any finite set / C IN the random measure 

ie{l,...,n}\/ 

converges weakly as n —> 00 to PPP(i^), the Poisson Point Process on 11 with intensity law v , for 
the usual vague topology on Radon measures. 

We refer to [531 Theorem 5.3, p. 138] for a proof of Lemma [2.51 Note that for i^j"' = a„/Uij 

it is a consequence of Lemma [2.41 (iii). In the case ^j"'' = an/Xij, where Xi,j is as in (jl.6p and 
(|1.7p . the above Lemma yields convergence to PPP{iya,e), where 

i^aAdx) = [M{,>o} + (1 - 0)t{,<o}] a\x\"-^dx , (2.11) 

2.5. Local weak convergence to PWIT. In the previous paragraph we have considered the 
convergence of the first row of the matrix a^^^X. Here we generalize this by characterizing the 
limiting local structure of the complete graph with marks an/Xij. Our argument is based on a 
technical generalization of an argument borrowed from Aldous [1 . This will lead us to Theorems 
[Q and below. 

Let Gn be the complete network on {1, . . . , n} whose mark on edge («, j) equals ^f^, for some 
collection {Q))i^i^j^n of i.i.d. random variables with values in K,, with ^"^ = Q]y We consider 
the rooted network (G„, 1) obtained by distinguishing the vertex labeled 1. 

We follow Aldous [U Section 3]. For every fixed realization of the marks (^[J), and for any 
B, F G M, such that {B"+^ - 1)/{B - 1) ^ n, we define a finite rooted subnetwork (G„, 1)^'^ of 
(G„, 1), whose vertex set coincides with a B-ary tree of depth H with root at 1. 
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To this end we partially index the vertices of (G„, 1) as elements in 

JB,H^Uto{h---,BYcMf, 

the indexing being given by an injective map (T„ from Jb.h to Vn {1, . . . , n}. The map a„ can 
be extended to a bijection from a subset of IN-'' to Vn- We set I0 = {1} and the index of the root 1 
is (T~^(l) = 0. The vertex v E Vn\l0 is given the index (fc) = a^^{v), 1 ^ fc ^ S, if ^"^^ has the 

fc**^ smallest absolute value among {£," j, j ^ 1}, the marks of edges emanating from the root 1. 
We break ties by using the lexicographic order. This defines the first generation. Now let Ii be the 
union of I0 and the B vertices that have been selected, li H ^ 2, we repeat the indexing procedure 
for the vertex indexed by (1) (the first child) on the set Vn\Ii. We obtain a new set {11, . . . , IB} of 
vertices sorted by their weights as before (for short notation, we concatenate the vector (1,1) into 
11). Then we define I2 as the union of Ii and this new collection. We repeat the procedure for (2) 
on Vn\l2 and obtain a new set {21, . . . , 2B}, and so on. When we have constructed {Bl, • • • , BB}, 
we have finished the second generation (depth 2) and we have indexed {B^ — 1)/{B — 1) vertices. 
The indexing procedure is then repeated until depth H so that {B^^^ — 1)/{B — 1) vertices are 
sorted. Call this set of vertices V^'^ = (^nJB.H- The subnetwork of Gn generated by V^'^ is 
denoted {Gn, 1)^'^ (it can be identified with the original network G„ where any edge e touching 
the complement of Vj^'^ is given a mark Xe — 00). In (Gn, l)^''^, the set {ul, • • • , u_B} is called 
the set of children or offsprings of the vertex u. Note that while the vertex set has been given a tree 
structure, (G„, 1)^'^ is still a complete network. The next proposition shows that it nevertheless 
converges to a tree (i.e. all circuits vanish, or equivalently, the extra marks diverge to 00) if the 
SJ^j satisfy a suitable scaling assumption. 

Let (T, 0) denote the infinite random rooted network with distribution PWIT(i/). We call 
(T, 0)^'^ the finite random network obtained by the sorting procedure described in the previous 
paragraph. Namely, (T, 0)^'^ consists of the sub-tree with vertices of the form u S Jb,h, with 
the marks inherited from the infinite tree. If an edge is not present in (T, 0)^'^, we assign to it 
the mark +00. 

We say that the sequence of random finite networks (G„, 1)^''^ converges in distribution (as 
n — 00) to the random finite network (T, 0)^'^ if the joint distributions of the marks converge 
weakly. To make this precise we have to add the points {±cx)} as possible values for each mark, 
and continuous functions on the space of marks have to be understood as functions such that the 
limit as any one of the marks diverges to +00 exists and coincides with the limit as the same mark 
diverges to —00. The next proposition generalizes [1] Section 3]. 

Proposition 2.6 (Local weak convergence to a tree). Let {£,i'j)i^i^j^n be a collection of i.i.d. 
random variables with values m R := RU {±cx)} and set ^"^ = S^^. Let u be a Radon measure on 
R with no mass at and assume that 

nP{^i2 e ■) as n^oo. (2.12) 

Let Gn be the complete network on {1, . . . , n} whose mark on edge {i,j) equals ^^j. Then, for all 
integers B, H , as n goes to infinity, in distribution, 

(G„,i)^-^^(r,0)^'^. 

Moreover, i/7i,72 are independent with common law PWIT(i^), then, in distribution, 

{{Gn, if'", {Gn,2f'") {{Tl,0f'", {%, 0f'"). 

The second statement is the convergence of the joint law of the finite networks, where (G„, 2)^'^ 
is obtained with the same procedure as for (G„, 1)^'^, by starting from the vertex 2 instead of 1. 
In particular, the second statement implies the first. 

This type of convergence is often referred to as local weak convergence, a notion introduced by 
Benjamini and Schramm [8], Aldous and Steele 0, see also Aldous and Lyons [2]. Let us give some 
examples of application of this proposition. Consider the case where = 1 with probability X/n 
and = 00 otherwise. The network G„ is an Erdos-Renyi random graph with parameter X/n. 
From the proposition, we retrieve the well-known fact that it locally converges to the tree of a Yule 
process of intensity A. If ^"^ = nYij, where l^.j is any non negative continuous random variable 
with density 1 at 0+, then the network converges to PWIT(i?i), where ii is the Lebesgue measure 
on [0, cxd). The relevant application for our purpose is given by the choice — {an/Xij), and 



12 



CHARLES BORDENAVE, PIETRO CAPUTO, AND DJALIL CHAFAI 



V = Vafi, where Xi,j are such that \Xij\ € and (jl.7p is satisfied, and Va^e is defined by (j2.11l) . 
Note that the proposition appHes to all a > in this setting. 

Proof of Provosition [KBl We order the elements of Jb.h in the lexicographic order, i.e. -< 1 -< 
2^---^i3^11^12^----<S...B. For v e Jbm, let Ov denote the set of offsprings of v 
in (G„, 1)-^'^. By construction, we have I0 = {1} and /v = cr„(Uw^vOw)- At every step of the 
indexing procedure, we sort the marks of the neighboring edges that have not been explored at an 
earher step {!,••• , n}\Ii, {1, . . . , n}\/2, Therefore, for all u, 

(C„(u),i)»^/u ^ 1^1(^11-11^1 ■ (2-13) 

Thus, from Lemma [2.51 and the independence of the variables we infer that the marks from a 
parent to its offsprings in (G„, 1)^^-'^ converge weakly to those in (T, 0)^'^. We now check that 
all other marks diverge to infinity. For v,w e Jb,h, we define 



'^v.w "^cr^ {v),o-„ (w) ■ 



Also, let {y^^, v,w g Jb,h} denote independent variables distributed as \£,i2\- Let E^'^ 
denote the set of edges {u, v} S Jb,h x Jb,h that do not belong to the finite tree (i.e. there 
is no k G {1, . . . , B} such that u = vfc or v = ufc). Lemma 12.71 below implies that the vector 
{|a;^ ,^1 , {v, w} e E^'"} stochastically dominates the vector 3^" := {y^.w ' {v, w} e E^'"}, i.e. 
there exists a coupling of the two vectors such that almost surely l^^v wl ^ 2/v wj fo'" {v,w} G 
E^'^ . Since Jb.h is finite (independent of n), 3^" contains a finite number of variables, and (j2.12p 
implies that the probability of the event {min^^ ^^^^b,h ^| ^ t} goes to as n — > 00, for any 
t > 0. Therefore it is now standard to obtain that if Xe denote the mark of edge e in T^'^, the 
finite collection of marks (a;")egj^ hxJb h converges in distribution to {xc)eejB h'xJb h as n — 00. 
In other words, (G„, 1)^'^ converges in distribution to (T, 0)^'^ . 

It remains to prove the second statement. It is an extension of the above argument. We consider 
the two subnetworks (G„, 1)^-^^ and (G„, 2)^-^ obtained from (G„, 1) and (G„, 2). This gives rise 
to two increasing sequences of sets of vertices /v,i and /v. 2 with v G Jb.h and two injective maps 
Cn,!, cr„^2 from Jb,h to {?1, • • • ,n}. We need to show that in distribution, 

((G„, 1)^^^, (G„, 2f^") ^ ((Ti, 0)^^^, (r2, 0f-"). (2.14) 

Let Vj^-^ = (Tn,i{JB,H) be the vertex set of (G„,i)^'^, i = 1,2. There are 

BH+1 _ 1 



G := 



B-1 



vertices in V^f^^ hence the exchangeability of the variables implies that 

P(2 e < ^. 

Let G„ = G„\F^{^, the subnetwork of G„ spanned by the vertex set V\V^i^ . Assuming that 
2(5^+1 _ - 1) < n and 2 ^ V^f , we may then define (G„, 2)^'^. If 2 e V^f , (G„, 2)^^" 

is defined arbitrarily. The above analysis shows that, in distribution, 

((G„, 1)^^^, (G„, 2f-") ^ ((Ti, 0)^^^, (r2, 0)^'"). 
Therefore in order to prove (j2.14p it is sufficient to prove that with probability tending to 1, 

K.Tnv;y = 0. 

Indeed, on the event {V^^{" n V^f = 0}, (G„,2)^'^ and (G„,2)^'^ are equal. For v e Jbm, 
let Ov,2 denote the set of offsprings of v in (G„, 2)^'^. We have 

/v,2 = {2}|JUw^vOw,2, 

and 

B...B 

nv,^f n vlf ^ 0) < P(2 e v,^f) + ^ P(0v,2 n v^f + \ v^f n 1^,2 = 0). 
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For any u, v e Jb,Hj Kfi^ ^ ^v,2 = 0, then (T„,2(v) is neither the ancestor of an,i{u), nor an 
offspring of cr„_i(u). From Lemma 12. 71 below we deduce that ^(u)^. ^(v) I gi'^^'^ ^fi^ '^-^v,2 = 
dominates stochastically |^"2|j and is independent ofthei.i.d. vector (|^^^^ ^^^-^ k\)ke{i n}\{v^-"ui 2 
with law |'^"2l- It follows that 

B 



P(a„,i(u)eOv,2lK:i n/v,2 = 0)<- „ I, 

n — o — \J-v.2 



Therefore, 



¥{0^,2 n v^f ^0\v^f n 1^.2 = 0) < ^PK,l(u)eOv.2lK^Tn-^v,2 = 0) 



CB 



Finally, 



nVnf n v^f ^ 0) ^ 



2C 

C C^B 



n 



- 2C" 



which converges to as n — >■ 00. □ 

We have used the following stochastic domination lemma. For any B, H and n let £^'^ denote 
the (random) set of edges of the complete graph on {1, . . . ,n}, such that cr~^(j')} is 

not an edge of the finite tree on Jb,h- construction, any loop belongs to E^'^- Also, for 
u 7^ on the finite tree, let ^(u) denote the parent of u. 

Lemma 2.7 (Stochastic domination). For any n E IN, and B,H E'iN such that 

BH+i _ 1 



B-l 

the random variables 

stochastically dominate i.i.d. random variables with the same law as law |C"2l- Moreover, for every 
u E Jb,h , the random variables 

{IC(u)J, *€{l,...,n}\a„(g(u))}, 
stochastically dominate i.i.d. random variables with the same law as law |^"2l- 



Proof of Lemma \2. 7[ The censoring process which deletes the edges that belong to the tree on 
Jb,h has the property that at each step the B lowest absolute values are deleted from some fresh 
(previously unexplored) subset of edge marks. Using this and the fact that the edge marks are 
i.i.d. we see that both claims in the lemma are implied by the following simple statement. 

Let Yi, . . . denote i.i.d. positive random variables. Suppose to = ni + ■■■ni, for some 
positive integers £, ni, . . . , n^, and partition the to variables in i blocks . . . , of rii, . . . , 
variables each. Fix some non negative integers kj such that kj ^ Uj and call ql, . . . ,qj.., the 
(random) indexes of the kj lowest values of the variables in the block P (so that F^i is the lowest 
of the Yi, . . . , Yn^ , y^i is the second lowest of the Yi, . . . , Y^ and so on). Consider the random 

index sets of the kj minimal values in the block, J-' := ^^LiWi}^ and set J = uj^^ J-'. If kj = 
we set = 0. Finally, let Y denote the vector {Yi , i — 1, . . . ,m ; i ^ J}. Then we claim that Y 
stochastically dominates m — X]j=i % i-i-'^- copies of Yi. 

Indeed, the coupling can be constructed as follows. We first extract a realization yi, . . . , of the 
whole vector. Given this we isolate the index sets J^, . . . , within each block. We then consider 
two vectors Z, V obtained as follows. The vector Zi = {zl, . . . , z^^^^,^^, zf, . . . , z^^^,^^, z^^^_^^) is 
obtained by extracting the ni — ki values z\, . . . , z^_^_f._^ uniformly at random (without replacement) 
from the values (in the block /^), the n2 — ^2 variables z^, . . . , in the same 

way from the values ym+i, . ■ ■ ,yni+n2 (in the block I^), and so on. On the other hand, the 
vector V = {vl, . . . ,vl^^_i^^,vf, . . . ,vf^^_f.^, . . . ,v^^_f.J is obtained as follows. For the first block 
we take vj, i — — fci equal to zj whenever an index i E \ was picked for the 

vector zl, . . . , z^^_f._^ and we assign the remaining values (if any) through an independent uniform 
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permutation of those variables yi, i \ which were not picked for the vector z;[, . . . , 2:^^_j,^. 
We repeat this procedure for all other blocks to assign all values of V. By construction, V ^ Z 
coordinate-wise. The conclusion follows from the observation that Z is distributed like a vector of 
m — X]j=i i-i-d- copies of Yi, while V is distributed like our vector Y. □ 

2.6. Proof of Theorem [2731 

Proof of Theorem \ 2. S]f i}. Let ly = v^fi be as in (|2.1ip and let {Ta,l2S) be a realization of the 
PWIT(j/). The mark on edge (v,vfc) in Ta is denoted by a;(v,vfc) or simply x^k- By definition, 
we have a;(v,w) = c» if v and w are at graph-distance different from 1. In particular, if we set 
Uv — sign(a;v)|a;v|", then the point sets Sv = {j/vfelfe^i are independent Poisson point processes 
of intensity £g = 61^x>a}dx -I- (1 — 6)t^j.^Qjdx. We may thus build a realization of the operator 
T on 7^, cf. (j2.3p . Let G„ be the complete network on {1, . . . ,n} whose mark on edge («, j) is 
:= Un/Xij. Next, we apply Proposition 12.61 For all B,H, (G„,l)'^'^ converges weakly to 
(7^,0)^'^. Let (T,f be the map (T„ associated to the network (G„, 1)^'-^ (see the construction 
given before Proposition 12. 6p . From the Skorokhod Representation Theorem we may assume that 
(Gn, 1)^'^ converges a.s. to {Ta, 0)^'^ for all B, H . Thus we may find sequences iJ„ tending 
to infinity, such that — 1)/(_B„ — 1) ^ n and such that for any pair u, v e IN-^ we have 

^(5 (u) a (v)) ~^ ^(m.v) a.s. as n ^ oo, where (t„ := cr^"'^". The map an can be extended to a 
bijection IN-^ — > M-^. It follows that a.s. 

(<5u,5-i(a-iX)?„5v) = ^ > ^— = (5u,T5v) . (2.15) 

^(S„(u),5„(v)) ^(u,v) 

Fix V S IN-^, and set V'n '-^ {o-n^ X)^,^^^ ■ To prove Theorem 2.3(i) it is sufficient to show that 
~^ T(5v in L-^iM^ ) almost surely as n — > oo, i.e. 

^((5u,C)-('5u,T<5v)f^O. (2.16) 

u 

Since from (|2.15p we know that (Juj'/'n) ~^ ((^ujTJv) for every u, the claim follows if we have 
(almost surely) uniform (in n) square-integrability of (((5u, '0n))u- This in turn follows from Lemma 
12.71 and Lemma l^^ i) . The proof of Theorem 12. 3f i) is complete. □ 

Proof of Theorem \ 2. S\f ii). We need the following two facts: 

lim — = 1 in probability, (2-17) 

n— )-oo ri 

and there exists S > such that 

liminf min — > 5 a.s. (2-18) 

n— >-oo l^i^n n 

Clearly, (j2.17p is a law of large numbers and holds actually a.s. (recall that for a > 1 we assume 
the mean of Uij to be 1). Let us establish the a.s. uniform bound (|2.18p . For every e > 0, there 
exists R> such that W,{Uijl^ij. j<R}) ^ 1 ~ e. If we define pf = ^i,j^{Ui j<R}^ then 

Pi 

lim inf min — ^ lim inf min — . 

n— >oo l^j^Ti n n— >oo l^i^n 71 

Therefore (|2.18|) is implied by the uniform law of large numbers in |lll Lemma 2.2], applied to the 
bounded variables Ui^jlfjj. .^j^. 

Next, we claim that for all u e W-'', in probability 

lim ^ 1. (2.19) 

n— ^oo Tl 

To prove this we first observe that by Lemma 12.71 and fj2.17l) we have in probabihty 

hmsup - ^^.(u).^.(.(u))) ^ ^ 

n— f oo ^ 

On the other hand C/5„(u),5„(g(u)) is stochastically dominated by the maximum of n i.i.d. variables 
with law Uij. The latter converges in distribution on the scale a„ (cf. Lemma l2.4f i)) and we know 



SPECTRUM OF MARKOV CHAINS ON COMPLETE GRAPH WITH HEAVY TAILS 



15 



that a„/n — >■ 0. It follows that in probability limsup„_j.oc Pan(-a)h^ ^ 1- Next, we can estimate 

'ie{l,...,n}\/u 

Now, observe that if u e IN^ belongs to generation /i, then the set /u contains at most 0{B^) 
elements, while n is at least of order B^" , where i3„,iJ„ are the sequences used in the proof of 
Theorem \2.3[ i). In particular, it follows that |/u| — o{n) and therefore (|2.13p and (j2.17p imply 
that liminf„_i.oo P5„(u)/"- ^ 1 in probability. This proves ()2.19p . 

Thanks to (j2.19l) . from the Slutsky lemma and the Skorokhod Representation Theorem, we may 
also assume that for each v e fi-^ , p-^^(y)/n converges a.s. to 1. We need to show that for each 
V € IN-^, (j2.16p holds with the new vector i/),^ := a~^{KnS)an5vj 

VP5„(v)P?„(w) 

Thanks to (j2.18p . ((<5w, V'^))w is uniformly square-integrable (cf. the proof of (|2.16l) '). and all we 
have to check is that ((5w,'0n) ^ (^w,T(5v))^ — ^ for fixed w. Here T is the operator appearing 
in the proof of Theorem 12. Sf i) above, now with the choice 9 = 1. We have 

{{5^,rn)-{K,T5^)f 

< 2 (a,T^C^5„(w),5,.(v) (1 - "/\/^5"(v)^r.(w)))^ + 2(a,7^[/5^(^)^5^(v) - {5^,T5^)f . 

The second term above converges to zero as in the proof of point (i). For the first term we use 
Pcf,i(v)/" — >■ 1 and P5„(w)/^ 1- This proves point (ii). □ 

Proof of Theorem \2.S\f iii). The setting is as in the proof of point (ii) above, but now a £ (0, 1). 
We build the operator S on the tree Ta as in (I2.6p . We need to prove that for any v e IN-^, a.s. 

^((<5w,O-('5w,S<5v))'^0, (2.20) 

W 

with '0^ := a~^SanSv, i-e. 

\/^S,i(v)P5„(w) 

Let us first show that for any v, w e IN-^ we have a.s. 

C^5?„(w),5„(v) ((5w,T(5v) 



^ . . (5w,S(5v). (2.21) 

V^S„(v)P5„(w) ^p(v)p(w) 

Multiplying and dividing by a„ and using (j2.15p with = 1, we see that (|2.2ip holds if 

«nV5„(v) ^p(v), (2.22) 
almost surely, for every v e fi-^ . In turn, (12.221) can be proved as follows. Let fc € M, and consider 



the tree with vertex set Jk,k, obtained as in Proposition 12.61 with B = H = k. Since J^^k is a finite 
set, for any v, (I2.15P implies that a.s. 

By Lemma [^771 and Lemma l^^ ii). X^u^/^ k '^"^ ^cf„(v),CT„(u) a-s. converges uniformly (in n) to as 
k goes to infinity. This proves ((^^ and 

Once we have (I2.2ip . to conclude the proof it is sufficient to show that a.s. 

lim sup (('5w,O)' = 0. (2.23) 

However, using and the simple bound ((^wjV'n))^ ^ , we have that again 

follows from an application of Lemma [^771 and Lcmma l^TW ii). This completes the proof of Theorem 
[OJiii). □ 
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2.7. Two-points local operator convergence. In the proof of the main theorems, we will need 
a stronger version of Theorem 12.31 Define the 2n x 2n matrices 

X ®X and S ® S, 

where "©" denotes the usual direct sum decomposition: X (B X{(j)i,(j>2) = (X(l)i, X(f)2), for n- 
dimensional vectors 01, 02- As for the limiting operators, we realize them on the Hilbert space 
L^iV) © L^iV) with V = K^. We consider two independent realizations Tj;, of the PWIT(^e), 
and call Ti, Si,T2, S2 the associated operators as in Section [231 We may then define 

Ti®T2 and Si ® S2. 

By Proposition EH {{Gn,l))^'" , {Gn,^)^'") converges weakly to {{T^ , 0)^^" , {T^ , 0)^^"). As 
before we can view the matrices X (B X and S' ® S' as bounded self-adjoint operators on L'^{V) ® 
L'^{V). Therefore, arguing as in the proof of Theorem 12.31 it follows that, in distribution, for all 
(01,02) e^?x 2?, 

a-^a-^X ® X(T„(0i, 02) ^ Ti ® T2(0i, 02) , 
where, tJ„ — cr^ ® ct^, and, as above, for i e {1,2}, cr^ is a bijection on IN-'', extension of the 
injective indexing map from to {1, . . . , n}, such that cr*i(0) = i. Analogous convergence results 
hold for the matrix S(D S. We can thus extend the statement of Theorem 12. 31 to the following local 
convergence of operators in L^{V) ® L'^{V). To avoid lengthy repetitions we omit the details of 
the proof. 

Theorem 2.8. As n goes to infinity, in distribution, 

(i) ifae (0, 2) then (a-^X ® a'^X, (1, 2)) ^ (Ti ® T2, (0, 0)), 

(ii) ifaG (1,2) and 9 = 1 then {unS ® KnS, {1,2)) ^ (Ti ® T2, (0, 0)), 

(iii) ifae (0,1) then {S®S, (1,2)) ^ (Si ® S2, (0, 0)). 

3. Convergence of the Empirical Spectral Distributions 

3.1. Markov matrix, a E (0,1): Proof of Theorem 11.41 Recall that S is a bounded self- 
adjoint operator on L'^{V), whose spectrum is contained in [—1, 1], cf. (|2.7p . The resolvents of S 
and S are the functions on C+ = {z S C : > 0}: 

i?(")(z) = (S*- z/)-i and R{z) = {S - zl)~\ 

For £ € IN, set 

pe:^{S0,S'S^) (3.1) 
Note that pe = {S0 ,^^60) p is the probability that the random walk on the PWIT associated 
to the stochastic operator K comes back to the root (where it started) after £ steps. In particular, 
Pe = for £ odd. We set pg = 1. Let ^0 denote the spectral measure of S associated to 60 (see 
e.g. [23l Chapter VII]). Equivalently, fj,0 is the spectral measure of K associated to the L^{V,p) 
normalized vector S0 :— 60/ ^ p(0), cf. (|2.7p . In particular, [I0 is a probability measure supported 
on [—1,1] and such that pg — ^\^x^ [i0(dx) , for every I. Since all odd moments vanish [i0 is 
symmetric. Moreover, for any z E C+ wc have 

{d0,R(z)d0)= / , 

J-l x~z 

i.e. {60 , R{z)60) is the Cauchy-Stieltjes transform of /i0. Recall that the Cauchy-Stieltjes trans- 
form of a probability measure ^ on R is the analytic function on C+ given by 

li{dx) 



rripiz) 



X — z 



The function characterizes the measure /i, |m^(z)| ^ {'^z) ^ , and weak convergence of /i„ to 
^ is equivalent to the convergence mp^{z) — > mp{z) for all z G C+. By construction 

itri?(")(z)= r^^=™,,(z), 
n J _i X — z 

where is the BSD of K, which coincides with the BSD of S. Using exchangeability and linearity, 
we get 

Ei?^"j'(z) = Eto^^(z) = TO]E^^.(z). 



SPECTRUM OF MARKOV CHAINS ON COMPLETE GRAPH WITH HEAVY TAILS 



17 



bmce 

i?'"'(^)i,i i^z)-'^ is bounded, we may apply Theorem 12.21 and Theorem 12.31 and obtain, 
for all z e C+, 

liui m-Ef^K{z) = ■m-Ef_,^{z). (3.2) 

n— f oo 

We define 

Next, we shall prove that, for all z G C+: 



We have 



lim E |mp^(z) - mE^0(2:)| = 0. 



E|mp^(z) - mE^^{z)\ < E |m^^(z) - Emp^(z)| + |m]Ep^(z) - mE^^{z)\ 



(3.3) 



On the right hand side, the second term converges to by 



The first term is equal to 



E 



1 " 



By exchangeability, we note that 




Ei?('\) 



n{n — 1) 



E 



-"-1,1 



Ei?("i' 



R 



in) 



Ei?, 



in) 



e!('1 - Ei?: 



(.n)\ ( T,{n) 
1,1 ) \ "-2,2 



Ei?, 



(n) 



n(9z)^ 

Theorem [22] and Theorem [2]8] imply that (i?i_i(2;), i?2,2(z)) are asymptotically independent. Since 
these variables are bounded, they are also asymptotically uncorrelated, and p.3p follows. 

Finally, observe that the sequence of measures ^xk is a.s. tight. Therefore the convergence p.3p 
is sufficient to establish a.s. convergence of ^.k to /Iq. This completes the proof of Theorem ll.4l 



3.2. I.i.d. matrix, a G (0,2): Proof of Theorem [TT2l Set A,, 

the Cauchy-Stieltjes transform 



"^X. For z e C+, we define 



"^A„(^) 



d/^A„ (x) 



k=l 



where 



i?(")(z) = (v4„-z/)-i, 

is the resolvent of A„. By exchangeability, WimA^{z) — Ei?J"^^(z). From Proposition IA.2I we know 
that T is self-adjoint. Therefore from Theorem 12.21 and Theorem 12.31 we infer 

Etoa„ {z) E/i(z) , h{z) ((50, (T - ziy^5<z) . (3.4) 

As in the proof of Theorem ll.4l we may write E/i(z) — Em^^, — rriEpj, , that is the Cauchy-Stieltjes 
transform of the expected value of the random spectral measure fiiz associated to T at the root 
vector (50. From (j3.4p we obtain the weak convergence of E/i^,, to /Xq, := E/i0. To obtain a.s. weak 
convergence of /iA„ to from Lemma IB. II it suffices to prove the convergence of Cauchy- 
Stieltjes transforms as in p.3p . This in turn is obtained by repeating word by word the argument 
in the proof of Theorem 11.41 

Thus, we have obtained ^.a^ fJ-a almost surely. Since the operator T only depends on the 
two parameters a and 0, where the latter is defined by (|1.7|) . the LSD fia might still depend on the 
parameter 0. However, the fact that fia is independent of 9 follows from Lemma 14.21 below, which 
implies in particular that the values m^ji^^[it) = E[ft,(i<)], t > 0, are uniquely determined by and 
therefore by analyticity, all values 771^^(2), z £ C-|- are uniquely determined by a. This ends the 
proof of Theorem 11.21 

We remark that in the proof of Theorem 11.21 one can avoid establishing (|3.3|) plus almost sure 
tightness (Lemma [BTTfi)) as we do above. Namely, the convergence of expected values E/z^-ij^ — 
is sufficient. This follows from an apriori concentration estimate; see [12j . However, we did 



18 



CHARLES BORDENAVE, PIETRO CAPUTO, AND DJALIL CHAFAI 



that piece of extra work here since we need it anyway in the case of Markov matrices, where the 
mentioned concentration estimate is not available. 

3.3. Markov matrix, a G (1, 2): Proof of Theorem 11.31 The proof given above for the matrix 
An = On^^ applies without modifications to the new matrix An '■— KnS, where Sij — . In 
particular, we use Theorem l2.3f M). Theorem \2.8( ii) and Lemma IB.IT m) to obtain the a.s. weak 
convergence of fiA„ to fia ~ E/i^, where fi0 is the random spectral measure of T at the root. This 
ends the proof of Theorem 11.31 

3.4. Markov matrix, a = 1: Proof of Theorem 11.51 Suppose now that a — 1 and set 
Wn — Jq"^ ^J-'idx) and k„ = na^^w„. A close inspection of the proof of Theorem l2.3r M) and 
Theorem 11.31 reveals that all arguments used for a E (1,2) can be applied to the case a = 1 
without modifications except for the two estimates (|2.17p and (|2.18p . which have to be replaced 
by p.Sp and p.6p below respectively. For p.6p we shall use the hypothesis (jl.lOp on w„. Let us 
start by proving that, in probability 

lim = 1. (3.5) 

n-»-oo nWn 

We recall that, for fixed i, a~^{pi — nWn) converges in distribution to a 1-stable law, see for instance 
[201 Theorem 1]. Therefore it suffices to show that k„ = a^^nwn — cx). To see this we may argue 
as follows. Observe that, for any e > 



1=1 i=l 



where Vi ^ V2 ^ • • • are the ranked values of J7ij, j — 1, . . . ,n. From LemmajHiH*) the right hand 
side above, for any e > 0, converges to ^J2i^i^{e^xi^i}i where the Xi are distributed according 
to the PPP with intensity x~'^dx on (0,oo). While this sum is finite for every e > it is easily 
seen to diverge (logarithmically) for e — > 0. This achieves the proof of p.Sp . 

Next, we claim that if w„ satisfies (jl.lOp . then there exists 5 > Q such that, a.s. 

liminf min ^' '^5. (3-6) 
To establish p.6p . let us define &„ — ayn^\ so that Wj{Ui,i\{u^ i^b„}) = wyn^\ and 

n 

Pi ^ "Sn := ^ ^l,»l{!7i..^b„}- 
i=l 

From the union bound. 



Ki 



mm 



l^i^n nWn ' \nw„ 

From Borel-Cantelli Lemma, it is thus sufficient to prove that for some 5 > 

nP (Sn < Snwn) < 00 (3.7) 

By assumption, there exists S > such that for all n large enough, w^ntj 2dwn- We define 

n 

Vi = C/i,il{c/i,.^&„} - win-l and Sn = ^ Vi. 

1=1 

Note that V,Vi = E5„ = 0. We get for all n large enough 

P {Sn < SnWn) — P {Sn < SnWn — nw^n^j) ^ P {S71 < —SnWn) ■ (3-8) 

By construction, Wn is slowly varying and a„ — L(n)n where L{n) is slowly varying. Hence 
\Vi\ ^ max(u'L„ej , 6„) — L{n)'nf where i(n) is another slowly varying sequence. By the Hoeffding 
inequality, we get from p.Sp 

P {^n < -SnWn) exp ( - ^^^^) = exp [~L{nW~'^) . 

where L{n) is a slowly varying sequence. Since e < 1/2 we obtain (13.71) and thus p.6p . 
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4. Properties of the Limiting Spectral Distributions 

Recall that /ic is characterized by the Cauchy-Stieltjes transform m^i^(z) = E/i(z), z e C+, 
where h,(z) is the random variable h^z) = (50, (T — 2/)~^(50), cf. p.4[) . The main novelty in our 
analysis of the LSD [i^ with respect to previous works [71 [S] is that we can work here with the 
distribution of hi^z^ rather than only with its expectation. 

4.1. Recursive distributional equation. The symbol = stands for equality in distribution. The 
following result is at the heart of our analysis of the LSD [1^. 

Theorem 4.1 (Recursive Distributional Equation). For all z G C+, the random variable 

h{z) = {d0,{T^~ziy^s^) 

satisfies to h{—z) = —h{z) and 

hiz)^-(z+J2M^)] , (4.1) 
\ fceK / 

where {hk)k<£w{z) are i.i.d. with the same law of h{z), and {(,k}keK is an independent Poisson 
point process with intensity ^x^^^^dx on {0,oo). 

Proof. Since the PWIT is bipartite, the property h{—z) — ~h{z) is a consequence of Lemma 
lA.ll We are left with the RDE (j4.ip . This can be interpreted as an operator version of Schur 
complement formula (see e.g. Proposition 2.1 in Klein |19j for a similar argument). Denote, as 
usual, by fc e IN the descendants of the root and let T^'^-' denote the subtree rooted at k (the set 
of vertices of T^'^-' is then kfi^). We have the direct sum decomposition IN^^ = {0}[JUkk¥i-^ ■ We 
define T^'^' as the projection of T on kK^ . Its skeleton is thus T^''\ Finally, define the operator 
U on I? by its matrix elements 

uk {S0,U6k) = (4,U(50) = ((50, T4) 

for all fc € ISf (offsprings of 0) and {6u, U(5v) = otherwise. In this way we have 

T = U + T with T = 0T('=^ 

fceK 

As T, each T'^'^^ can be extended to a self-adjoint operator, which we denote again by T'^'^^. 
Therefore T is self-adjoint. We shah write R{z) = (T - zl)'^ and R{z) = (T - zl)'^ for the 
associated resolvents, z E C+. These operators satisfy the resolvent identity 

R{z){T -T)R{z) ^ R{z) - R{z) . (4.2) 

Set Ru,viz) :— {Su, Riz)Sv) and i?u,v(-z) :— {Su, R{z)Sv). Observe that R0,0{z) = —z~^ and that 
the direct sum decomposition IN^^ = {0} IJ UkkK-^ implies Rk^i{z) = for k ^l. Similarly we have 
that R0^k{z) = = Rk,0{z) for every fc e E^. From (j4.2p we then obtain, for A; S M: 

Rk,k{z)ukR0,0{z) = ~Rk.,0{z) ■ 

It follows that 

((50, i?(z)(T - T)R{z)50) = R0,0{z>kRk,0{z) = - ^ R0,0{z)RkMz)ulR0,0{z) . 

feeiN ken 

From (|4.2p we then conclude that 



R0.0{z) 



R0,0{z) 



1 - R0,0{z) 'Eken Rk,k{z)ul 
Or, using R0^0{z) — —z~^: 



R0,0(z) 



z + ^ Rk,k{z)ul J 
\ ken / 



Then (|4.ip follows from the recursive construction of the PWIT: T^^'^ are i.i.d. with distribution 
T and therefore Rk,k{z) are i.i.d. with the same law of R0^0{z), for every z G C-|-. □ 
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Concerning the uniqueness of the sokition to the RDE (j4.ip we can establish the foUowing useful 
result. For z — it, with t > 0, the identity, h{—z) = —h{z) reads ^h{it) = 0. Thus, the equation 
satisfied by g{it) = ^h{it) ^ is 

3(zt)= U+ Va5fc(*i) I • (4.3) 



\ keTN / 



Lemma 4.2 (Uniqueness of solution for the RDE). For each t > 0, there exists a unique probability 
measure i** on 11+ , solution of 



Proof of Lemma \4.'2\ Set /3 = a/2. If {Yk) is an i.i.d. sequence of non negative random variables, 
independent of {C/tjfceiN, such that E[y]'^] < oo then it is well known that 

k k 

(see for example [27l Lemma 6.5.1] or (j4.5p below). This implies the unicity for Equation (j4.3p 
provided that the equation satisfied by W,[g{it)^] has a unique solution. Recall the formulas of 
Laplace transforms, for y ^ 0, rj > and < r/ < 1 respectively, 

y^'^ = T{ij)-^ x'^-^e^'^ydx and = r(l - 77)"^ / a;~''"^(l - e~^^)dT. (4.4) 

JQ Jo 

From the exponential formula we deduce that, with s ^ 0, 

Eexp (^-s^^kYk^ = cxp (l]J^ (e-^"^i - Ijpx'^'^dx 

= exp (-r(l - /3)s'^E[Yf]^ . (4.5) 
From Equation (|4.3p . 'E[g{it)^] is solution of the equation in y: 

V = — — / x^-^e-^'^e^^'^^'^^'^^ydx 

The last equation has a unique solution for any i ^ 0. Indeed, the function from IR,+ to IR+ 

1 



m Jo 

tends to as y — >■ oo and it is decreasing since 



r(/3) 

Thus (/3 has a unique fixed point. □ 



Before going into the proof of Theorem ll.61 we introduce some notation. Let P ^ a/2 as above 
and let fCa denote the set of probability measures on (0,oo) with finite (5 moment. We define the 
map ^ on probability measures on R_|_ U {oo}, where ^'(Q) is the law of 

^ = f E ^^^fe) ' (4.6) 
VfceiN / 

with {Yk, G IN) i.i.d. with law Q independent of S = {CfelfeeiN a Poisson point process on R_|_ of 
intensity j3x~^~^dx. 

Lemma 4.3. 4* satisfies the following 

(i) is a map from K-a to K-a- Let (Pn)neN and P in K-a, i/lim„_>.oo / x^dPn — J x^dP then 
*(P„) converges weakly to ^'(P) and lim„_>oo / a;^(i^'(P„) = / a;^d*(P). 

(ii) The unique fixed point of in ICa is the law ofl/S where S is the one-sided j3-stable law 

with Laplace transform Eexp(-tS') = cxp (^-i'^^r(l + /3)/r(l - /3)^, t > 0. 

(iii) E5-^ = (r(/3 + i)r(i - /?)))-i/^ 
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Proof of Lemma \4.3\ As in the proof of Lemma 14.21 we get 



E 



1 



m Jo 
1 



m Jo 

" pmJo 

= (r(/3 + i)r(i-/?)Er/)-i, 

(in the last hne we have used the identity zT{z) — T(z + 1)). Therefore, \l/ is a map from /Cq to 
JCa- Also as a consequence of (|4.5p : 

Eexp(-tZ-i) = exp(-i''r(l - /3)EFj^). 

Statement (i) follows from the continuity of the map x 1/a; in (0, oo). If Z is a fixed point of ^ 
then from the computation above EZ'' = (r(/3 + l)r(l - (3)))-^/^. Finally, from (g^]) we obtain 
for all t ^ 0, 



Eexp(-tZ-i) = exp(-i'5r(l - /3)EZ^) - exp ( -i^ J£ii±|| 



□ 



4.2. Proof of Theorem From Theorem 14. 1[ for z G C+, 

where h solves RDE (14. ip . Set f{z) = 3fi/i(z) and g{z) = For z = m + w G C+, / and g 

satisfy the RDE 

U + J2kikfkiz) 



and 



iu + EkUkiz)f + {v + j:^^k9kiz)f 



V + J2k^kgk{z) 



{u + Ek^kfk{z)f + {v + j:,ik9k{z)f' 



By construction, ^ g{z) ^ thus the law of g{z) is in ICa- If the stochastic domination of P 
by Q is denoted by P Q, we have 

/ \ -1 / \ -1 



g{z) K + X! ^kgk{z) ^'st X! ^kgkiz) 



(4.7) 



(In fact, we also have |ft-(z)| (X]fe ^kgk{z)) ^). Using the computation in Lemma 14.31 we obtain 
E5(z)^ (r(/3 + l)r(l - (3)W.g{z)^)-\ Thus 

E5(z)'5 «C ^ (4.8) 

Vr(/? + i)r(i-/3)) 

Again, the formula y^^ — r(77)^^ x'^^^e~^'^dx, for y ^ 0, 77 > 0, gives 



E 



1 

m Jo 



^^-i^-.^ni-mgi^fdx. 



(4.9) 



We now study the weak limit of g{u + iv) when v i 0, u G M,. Equation (|4.8p implies tightness, 
so let g{u + iO) be a weak limit. If this limit is non zero then 'Eg^{u + iO) > 0, and Equations 
(P?71) - (|i^ imply for all 77 > and u G R, 

limsup W,g^{u + iv) < 00. 

u-\-iv:v].0 
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Since Wih{z) is the Cauchy-Stieltjes transform of /ic, taking rj = 1, we deduce that fia is absolutely 
continuous with bounded density, see for example [ISl Theorem 11.6]. □ 

4.3. Proof of Theorem ILBI^ m). In view of Theorem 11.6], it is sufficient to show that 

Hnr%(^^)-r(l + i)(I|l±|)". (4.10) 

As above, (|4.8p implies the tightness of {g{it), t > 0). So let g{iO) be a weak limit. It is in ICa and, 
by continuity, g{iO) is solution of the RDE 

By Lemma Ol g{iO) = 1/S, and gives 



Eg(iO) = / e V r(i+« ^ _r , , , ^ 
Using the identity zr(z) = r(z + 1), we get (g^TUl). □ 

4.4. Proof of Theorem \1.6( iii). We start with a Tauberian-type theorem for the Cauchy- 
Stieltjes transform of symmetric probability measures. As usual, let denote the Cauchy-Stieltjes 
transform of a symmetric probability measure fi on R. Then, for all t > 0, rrif^^it) G iR+ and 

f°° t f°° t 

Lemma 4.4 (Tauberian-like lemma). If L is slowly varying and < a < 2, the following are 
equivalent: as t goes to +oo 

M((t,(XD))^L(i)i-" (4.11) 
QrUf^iit) - A{a)L{t)t~°'-^ (4.12) 

with A(a) = 2a j^^dx. 

Sketch of Proof of Lemma \4-4\ The proof is an adaptation of the proof of the Karamata Tauberian 
Theorem in |10[ pages 37-38]. Let Ai denote the set of symmetric measures on R such that 
/q min(l, a;^)/Lt(dx) < +oo. On M, define the transform 

f°° 2x^ 



Note that Sii{t) = 1 — iGm^(ii) = 1 + itm^{it). Recall that the Cauchy-Stieltjes transform 
characterizes the measure. Thus if for all < > 0, (5/i,i(t))„g]N converges to 5/i then (/i„)„g]N 
converges to ^ over all bounded continuous function with outside the support. Now, assume that 
(|4.12p holds, namely 

Sn{t) - A{a)L{t)t'°'. (4.13) 
Since lim2;__j.oo L{tx) / L{t) — 1, we deduce that for alH > 0, as 2: — > 00 

Sii{xt) 



L{x)x~°' 



A{a)t-°'. 



The left hand side is the S transform of the measure fJ-xidy) — fi{xdy)/{L{x)x ") while the right 
hand side is the S transform of /ioo(rfj/) = a\y\~°'~^dy, thus 

r( N -a = Aix((l,oo)) ^ /ioo(l,oo) = 1. 

lj\X JX 

We get precisely (|4.1ip . The reciprocal implication can be proved similarly, see [TUJ pages 37-38] 
(it is straightforward for L{t) — c, the case that we will actually use). □ 
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We now come back to the RDE gj]) and define Q{t) = E[g(ii)^]. From (|43)) . we have a.s. 
tg{it) ^ 1. Note also, from a.s. J2kikgk{'it) < t~^J2k^k, that a.s. \imt^+oDtg{it) = 1. The 
dominated convergence Theorem leads to 

lim t^^Qit) = 1. 



(4.14) 



Moreover, as already pointed in Lemma 14.21 



We deduce, with C(t) = {tQ{t)^/'^y^/^ , that 



E 



t' + tQity/i'Ek^k 
""^'^^{tcitW + Ekik 



(4.15) 



where C{Y) is the law of 



and e is independent of {^k}k, P(e 




1/2. We have 



By g3]), as s ; 0, E exp(-s Cfc) = exp(-s'3r(l - /3)) - 1 - sf^T{l - P). Using [TUl Corollary 
8.7.1], we obtain V{Y,k£.k > t) ^ and 

-i^^A(a). Thus by 



By Lemma Ol 3m£(y)(zt) - 

3m^t^ (it) — i^"' 

Theorem II. 6f iM) now follows from Lemma 14.41 



-A(a). 



□ 



Remark 4.5. /n the proof of Lemma \4.S\ we have seen that the distribution of g{it) — Qh(it) was 
function of Q{i) — 'E[g^{it)] which satisfies the equation 



Qit) 



fpit,Q{t)). 



We could push further our investigation at t — and compute the derivative of Q at t = 0: 
Q'(0) = -fp+iiO, Q(0)) - r(l - /3)/2,3(0, g(0))Q'(0), with Q(0) = (r(/3 + l)r(l - /3)))-i/2. There 
should be no obstacle for computing by recursion the successive derivatives of Q(t) att~0. We 
would then obtain a series expansion of the partition function fj,a{{—oo,t)) in a neighborhood o/O. 

4.5. Proof of Theorem II. 7t Jlai ol £ (0,1). As in p.ip . let pe denote the return probability 
after £ steps starting from the root 0, for the random walk on the PWIT with transition kernel K 
given by (12. 5p . In particular, -fi — Ep£ is the l-th moment of the LSD Jla- 

Proof of Theorem \1.7\ (i). For the first part, we shall show that there exists 6 > such that for 
any e e (0, 1/2] and any n: 

j2n^Se^{l-e)^". (4.16) 

Theorem 11.71 (i) follows by choosing e — l/2n. To prove (j4.16l) we use the simple bound P2n ^ 
(K(0, 1)K(1, 0))", which states that to come back to the root in 2n steps the walk can move to 
the child with the highest weight, with probability K(0, 1), go back to the root, with probability 
K(l, 0), and repeat this n times. Taking expectation, it follows that 

72„ ;^E[(K(0,1)K(1,0))"] . (4.17) 
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Therefore (j4.16l) holds if the event 

A, = {K(0, 1)^(1- e) and K(l, 0) ^ (1 - e)} 

has probabihty at least 5e", for some 5 > Q and for any e € (0, 1/2]. 

Let (a;i)i denote the realization of the PPP at the root 0, i.e. xi > 2:2 > ■ • ■ are the points of 
a PPP on (0,00) with intensity measure ax^'^~^dx. We set </> := YlTLi ^"^^ I'^t 4/ denote an 
independent copy of (j). We can use the representation K(0, 1) = xi/(t) and K(l, 0) = xi/{xi +</>')• 
Therefore, 

P(A,) = P (.Ti ^ (1 - £)0 , a;i ^ (1 - e)(xi + 0')) 

= p(a;i^ (!-£)</>, '/''^ 

^ P (xi ^ (1 - , xi > e"^ , 0' s$ 1) . 

Let 5i := P((/) ^ 1) = Jq f{t) dt > 0, where /(t) denotes the density of 0. The function f{t) can be 
obtained from its Laplace transform, which is given by the known identity E[e^'"^] = g-r(i-Q)ii ^ 
u > (see ^22j Proposition 10], or (14.51) with (3 replaced by a and Yk = 1). Since cj)' is independent 
of (xi) we obtain 

P(^,) ^ 5i P {xi ^ (1 - e)(/. , xi ^ e-i) . 

To estimate the last quantity we observe that if a; is a size-biased pick from (xi) then xi ^ x. We 
recall that x is a random variable such that, given the sequence (xi) the probability that x equals 
Xi is Xi/(j). It is not hard to check (see e.g. [21] Lemma 2.2]) that the random variable x has a 
probability density on (0, 00) given by 



fit) 



X + 1 



dt, 



where f{t) is the density of the variable cj). Therefore, 



P {xi ^ (1 - e)(l) , xi > e^i) ^ P (5 > (1 - £)0 , x ^ 



dtf{t) 



(4.18) 



di/(i) 



(5i (1 -e) 



In conclusion, P(A£) > (1 - e) > i 5l e", and the claim (|4T6)) follows. 

It remains to show that liminfa/ii 72 > 0. If {xi), x, and (f> are as above and if cj)' is independent 
of the sequence (x^ ) and identical in law to the random variable (f) then 



72 =IE 



E 



4> Xi + 4>' 



E 



X + I 



ax 



l-a 



' fit) 

X + t 



dt dx 



Now, from the Laplace transform E[e ' 

fit) 



X + t 



dt 



r(i have the identity 



This gives 



/"OO /"OO 

72 = ar(2 -a) / 
Jo Jo 

aT{2 - a) 



r(i-a) 



00 nOO 



it 



ds dt . 



Finally, the desired result follows from the bounds (for absolute constants Ci , C2 > 0) 



00 /"OO 



[t + s)-2+" ds dt ^ e-M / (t + s)-2+" dsdi ^5 



1 — a 
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and 

/•oo /•! /"OO 

r(l - a) = / t^"e~* < / t-" dt+ e-* dt — 
Jo Jo Ji 1 ^ 



a 

□ 



Proof of Theorem \1.7\ (ii). It is convenient to make here the dependence over a explicit in ah the 
notations. In particular, for every a e (0, 1), we denote by Sq the operator S given by (12.61) . These 
operators are defined on a common probability space, and are self-adjoint in L'^(V). Moreover, 
it follows from Subsection 13.11 that fia — E/^q,0, where fia^0 is the spectral measure of Sq. at the 
vector 60 . By the dominated convergence Theorem, in order to prove that a 1— >■ /Xq, is continuous 
in (0, 1), it is sufficient to show that a.s. a 1— ?> fj,a,0 is continuous. From [221 Theorem VIII. 25(a)], 
it is in turn sufficient to prove that for all v £ a i-> SaS^ is a continuous map from (0, 1) to 
L^(F). From (j2.6p . for all u G V, the map a H> Sq,(u,v) is continuous. It thus remains to check 
the uniform square integrability of (Sq,(v, u))^gy. We start with the upper bound 

— 1/q — 1/a — 1/a 

Pa(v)pa(vfc) p«(v) 

Then, notice that for all a G (0, 1 — e), one has y^l^°' ^ max(l, y.^,^!^'^ '^''), and Pa(v) ^ 
min(l, y.^,/^'^ ^^). We may conclude by recalling that a.s. lim^ yvk/k = 1 and y^i > 0. □ 

Proof of Theorem ] 1.7\ (in). As in the proof of Theorem [TjT] (ii), we make here the dependence over 
a explicit in all the notations. It follows from Subsection 13. II 



where the expectation is over the randomness of the PWIT. We introduce for v G y, 

K.(v) 



-1/a -1/a 

2/vi 



\Z^fc>li/vfc Z^k^l yvk / 

By construction Kt(v) is a PD(a,0) random variable. Thus, by [211 Corollary 18], as a | 0, Vq(v) 
converge weakly to the deterministic vector (1, 0, • • • ). We may thus write, 

-1/a 



Vi +2/11 (l + ^a) 

where as a goes to 0, goes in probability to 0. We define U — l{yii>j,i} , so that J7 is a symmetric 
Bernoulh i.e. P([/ = 0) = F{U = 1) = 1/2. We have proved that in probability. 



In particular. 

Since ^J,a,0 is symmetric. 



liniKa(0,l) = l and \iuiKJl,0) = U. 

aiO aiO 



lim / x^^ iXn, 0(dx) = U. 



lim^Q,0 = ^S^i + (1 - U)So + ^Si. 

alO A A 



Taking expectation, we obtain the claimed statement on fia. □ 

5. Invariant Measure: Proof of Theorem 11.81 

We start with a lemma. Let {Xi, . . . , X„), Xi ^ ■ ■ ■ ^ Xn, denote the ranked values of pi, . . . , pn 
and recall the notion of convergence in the space A, cf. Section We use the notation 6„ :— am„ , 
where m„ = n(n + l)/2. 

Lemma 5.1. For any a G (0,2), the sequence b~^{Xi, X2, ■ ■ ■) converges in distribution to 
{xi,Xi,X2,X2, ■ ■ ■), where Xi > X2 > ■ ■ ■ denote the ranked points of the Poisson point process 
on (0,00) with intensity ax~°'^^dx. 
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Proof of Lemma \5A[ There are m„ = n{n + l)/2 edges, including self-loops. Let us denote by Ug 
the weight of edge e S {!,..., nin}. The row sums are given Pi = J2e- e3i^e- We write 0„ for 
the set of off-diagonal edges e, i.e. edges of the form e = {«, j} with i ^ j. Let Ue^ ^ ^ ■ ■ ■ 
denote the ranked values of the i.i.d. random vector ([/e)eeo„- Since there are m„ — n edges in 
On , an application of Lemma 12.4^ 1) yields convergence in distribution 

b-\Ue„Ue„...) ^ (X1,X2,...). (5.1) 

n— >oo 

Each Ci = {ui,Vi} e On identifies two row sums p„. and py.. Set = maxjpu^ — Ue^yPvi — Ue^}- 
Then, for every fc g IN and e > 0: 

lim P ( max > e fe„ ) = . (5.2) 

To prove this we use an estimate due to Soshnikov _261 . Let Bn denote the event that there exists 
no i G {1, . . . , n} such that 

{pi > bn and pi - maxC/ij- > bn }. 
j 

Then, from and [H Lemma 3], one has 

lim V{Bn) 1 . (5.3) 

Clearly, on the event Bn, if maxi^^^^ Ag ^ e bn, then Ue,, ^ bn ^ which has vanishing probability 
in the limit by ((5?T|) . This proves ((5T2|) . 

For simplicity, we introduce the notation i?2£-i = rnax{p„^, i?2£ = mm.{pu(., Pve}- There- 
fore (|5.2p and (jS.ip prove that 

bn^{Rl,R2,R3,R4---) {xi,Xi,X2,X2,...) . (5.4) 

n— )-oo 

It remains to show that for every fixed k: 

lim P (Ui^,^2fe{i?, ^ X,}) = . (5.5) 

n— ^oo 

By construction, we have Xi ^ Ri for i — 1,2. On the event Bn described above, to have 

Xi > Ri or X2 > i?2 implies that there exists an edge e 7^ ei such that Ue ^ C/ei — bn ^ ■ 
However, this event has vanishing probability by (|5.1I) and the fact that b^^ maxiUi^i — >■ in 
probability for all sufficiently small i5 > (indeed by Lemma [2.41 On^ maxiUi^i converges weakly 
to the Frechet distribution, see first comment after Lemma [2.4^ . Thanks to (|5.3|1 this shows that 
P(Xi > i?i or X2 > R2) ^ 0. Recursively, the probability of X2i+i > i?2i+i or ^2^+2 > i?2i+2 
on the event Bn H {Xj — i?^ , Vj = 1, . . . , 2i} vanishes as n — > 00. Indeed, at each step we have 
removed a row and a column corresponding to the largest off-diagonal weight and we may repeat 
the same reasoning as above. This proves (j5.5p as required. □ 

Proof of Theorem ] 1.8\f ii). Let us define m„ = n{n + l)/2. Observe that 

n n 

^ Pi = 2Sn + Dn where Sn := ^ Ue and Dn := ^ J7i,i . (5.6) 

i=l eeO„ 2=1 

Here, as in the previous proof 0„ denotes the set of off-diagonal edges. For a G (1,2), we have by 
the weak law of large numbers S'„/to„ — > 1 and Dn/n — )• 1 in probability. Therefore 

1 " 

lim > Pi — 2 , in probability. (5.7) 

n->-oo m„ ^ — ' 
i—l 

Theorem ll.Sr ii') thus follows directly from Lemma 15.11 and (j5.7l) . The same reasoning applies in 
the case a = I replacing the law of large numbers by the statement p.5|) which now gives (j5.7p 
with m„ replaced by m„z«m„ . D 

Proof of Theorem ] 1.8\f i). If C/ei ^ ^ ••• are the ranked values of the i.i.d. random vector 
(C/e)e60„ and Sn is their sum as in (j5.6p then by Lemma l2.4f ii). replacing n with m„, we have 

' Uei Ue2 \ d ( Xi X2 \ , . 
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where Xi > X2 > ■ ■ ■ denote the ranked points of the Poisson point process on (0, oo) with intensity 

Write Xi,X2, ... for the ranked values of row sums as in Lemma lSTTl so that pi = Xi/{2Sn+Dn), 
where I?„ , 5„ are as in (j5.6l) . Let 

Y _ ^21-1 Uee Y — ^2^ Uee 

Thanks to (|5.8p it is sufficient to prove that P(maxi^i^2fc > e) — , as n — > oo, for any fixed 
e > and fc e IN. This follows from the argument used in the proof of (j5.2p and (j5.5p . □ 

Appendix A. Self-adjoint operators on PWIT 

The following classical lemma was used in Section |3l If S is a self-adjoint operator on D{S) C 
L'^{V) with V countable, the skeleton of S is the graph on V obtained by putting an edge between 
two vertices (v, w) iff {Sv, S(5w) 7^ 0. 

Lemma A.l (Resolvent of self-adjoint operators on bipartite graphs). Let S be a self-adjoint 
operator on D{S) C L'^{V) with V countable. If the skeleton is a bipartite graph then for v ^ V , 
h{z) — {Sv, (S — zI)^^Sv) satisfies for all z € C+, h{—z) = -~h{z). 

Proof of Lemma \A.l\ Assume first that S is bounded: for all ^ V , ||S(5w|| =^ C. For \z\ > C, 
the series expansion of the resolvent gives 

However since the skeleton is a bipartite graph, all cycles have an even length, and for £ odd 
(JvjS^Jv) = 0. We deduce that for \z\ > C, h{—z) = —h{z). We may then extend to C+ this last 
identity by analyticity. 

If S is not bounded, then S is limit of a sequence of bounded operators and we conclude by 
invoking Theorem VIII.25(a) in P]. □ 

The arguments of Section |3] were crucially based on the following fact. 

Proposition A. 2. The operator T defined by (|2.3p is essentially self-adjoint. 

To prove the proposition, we start with a deterministic lemma. Let V = IN-'' denote the vertex 
set of the PWIT and let T) be the space of finitely supported vectors. We write u ^ v if u = vfc 
or V = ufc for some fc G IN^ (i.e. if u,v are neighbors) and u 7^ v otherwise. Let yl : 2? — > L'^{V) 
denote the symmetric linear operator defined by 

((5v, A(5w) = Wu.v = Wv,u , (A.l) 
and such that w^.v — whenever u 7^ v. 

Lemma A. 3 (Criterion of self-adjointness) . Suppose that there exists a constant k > and a 
sequence of connected finite subsets (S'„)„^i in V, such that Sn C Sn+i, ^nSn = V, and for every 
n and v Cz Sn, 

Then the operator A defined by (jA.l[) is essentially self-adjoint. 
Proof. It is sufficient to check that the only function ip £ D{A*) C L'^{V) such that 

A*(p — ±iip 

is (fi = (see e.g. [531 Theorem VIII.3]). A similar argument is used in P^, Proposition 3]. We 
deal with the case A*ip = iip, i.e. for all u e 

i(p{u) = ^ Wu,v</5(v). 
v~u 

Here we use the notation (p{u) — {6u, f). Taking conjugate, we also have for all u S F 

-i^(u) = ^ Wu,v^(v) = ^ i(;v,u^(v). 
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For any finite set S C V, wc deduce 

i E I^WI' = E ^(v)(A»(v) = E ^(v) E ^v,u^(u) 
ves ves ves u~v 

= E "^("^ E ^v,u^(v) + E ^(v) E "'v.u'^^cu) - E E ^v,u^(v) 

ueS v~u ves u~v: u^S uSS v~u:v^S 

= E + E E "'v,u<^(u) - E vi'^) E ^v,u^(v). 

ueS veS u~v: u^S u6S v~u: v^S 

We obtain a Green formula: 

E = E '^^"^^ E ^v,u<p(u) - E 'fi'y^) E w'v.u^cu), 

vGS veS u~v: U0S veS u'-v: u^S 

From Cauchy-Schwarz inequality: 

Ei'^wi' < Ei'^(^)i E i^'v.uibiu)! 

^ (Ei^wi')'^'(E( E K,uik(u)i)y^'. 

veS veS u~v: u^S 

Now take S = Sn- Prom the assumption of the lemma, using again Cauchy-Schwarz inequality, 

( E kv,uii<^(u)i)'</« E i^wi'- 

Since Sn is connected and the graph is a tree, if u ^ S'n and u ~ v then for any v' € Sn\^, then 
u 7^ v'. It follows 

E i^(v)p ^ v^( E E i^wi')'^'- 

ves„ ves„ ues= 

Therefore, 

E i^(v)r<KEk(v)i'- 

Since lim„ Sn = 1^, as n grows, the right hand side goes to 0, while the left hand side goes to \\ip\\2- 
We obtain ip = 0. □ 

Next, we need a technical lemma. 

Lemma A. 4. Let k>0, 0<a<2 and let < xi < X2 < ■ ■ ■ be a Poisson process of intensity 1 
on 1R+. Define = inf{t e IN : x^^^" < k}. Then Er^ is finite and goes to as k goes 

to infinity. 

Proof. First of all, the fact that is a.s. finite follows from the a.s. summability of Y^^=i 
We deduce also that a.s. there exists k > such that = 0. From monotone convergence, it 
remains to check that Er^ < oo. Let n ^ 1 and Sn = X^^Li ^fe ^^"^{xfe^n}- Prom Levy-Khinchin 
formula, ioi 9 > 0, 

Eexp(6'S'„) = exp (^J (e^^^'^" - l)dx 
As n goes to infinity, if ^ = o(n^/"). 

Hence, taking 9 = (2/a — l)n^/"~^, we deduce from Chernov bound, that for any integer n > no, 

¥{Sn > k)^ e-*'*Eexp(6'5'„) ^ Se"^"''""', 

where no > 1 and c = {'2/a — !)«;. Also recall (from Chernov bound) that if A'' is a Poisson random 
variable with mean n then for alH > 0, 

PfiV > i) < exp ( -Hog — 
V ne 
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Now if the event {r^ > t} holds then either the number of points of the Poisson process {xk)k^i 
in [0, n] is larger than t or Sn > We get for any integer n ^ no, 

We conclude by taking n — max(no, t/(2e)). □ 



Proof of Provosition \A.'A We apply Lemma fA.BI with A given by T, the operator defined by (|2.3p . 
For K > and v e IN-^, we define the integer 

oo 

T«(v) =inf{t ^0: J2 lyvfcr'/"<4- 
fc=t+i 

The variables (tk(v))v are iid and by Lemma [A. 4[ there exists k > such that Etk(v) < 1. We fix 
such K. Next, we give a green color to all vertices v such that Tk(v) ^ 1 and a red color otherwise. 
We consider an exploration procedure starting from the root which stops at red vertices and goes 
on at green vertices. More formally, define the sub-forest of the PWIT where we put an edge 
between green vertices v and vfc iff 1 ^ fc ^ Tk (v) . 

The sets Sn appearing in Lemma lA.3l are defined as follows. If the root is red, we set Si = {0}. 
If the root is green, we consider T^, the maximal subtree of that contains the root. It is a 
Galton- Watson tree with offspring distribution t„(v). Thanks to our choice of k, is almost 
surely finite. Let denote the set of vertices of , and consider the set of the leaves of 
. Note that is the set of vertices v € V^J such that for all 1 ^ fc ^ ''"k(v), vfc is red. Thus, 
when the root is green, we set 5*1 = Uvei^ {vfc : 1 < fc ^ Tk(v)}. By construction, the set Si 
satisfies the condition of Lemma IA.31 

Next, define the outer boundary of the root as {0} as d{0} = {1, • • • ,Tk(0)}, and for v 7^ 0, 
V = (ii, • • • ik), set 

^{v} = {{h,--- ,ifc-i,ifc + 1)} U • • • ,ik,l),--- - Jk^Ti^i-v))} ■ 

For a finite connected set S', its outer boundary is defined by 

95=(U9{v})\5. 

ves 

To define the set S2, suppose that dSi = {ui,--- ,u„}. The above procedure defining for 
the PWIT rooted at can be now repeated for the subtrees rooted at wi, . . . , u„ to obtain sets 
5*1 (ui), . . . , 5*1 (u„). We can then define 5*2 = S'l IJ Ui^i^nSiiui). Iterating this procedure, we may 
thus almost surely define an increasing connected sequence (Sn) of vertices with the properties 
required in Lemma I A. 31 □ 

Appendix B. Tightness estimates 

Let X and K be the matrices defined by (|1.6p and (jl.ll) respectively. Recall that, when a ^ 1 
we set Kn — nwna~^^ where Wn = 1 if a > 1 and Wn — /q " xC(dx) if a = 1. 

Lemma B.l. 

(i) For every a € (0,2), the sequence fJ-^-^x '"^ i^Q^t- 

(ii) For every a € [1,2), the sequence fJ.K„K is a.s. tight. 

We first recall a classical lemma on truncated moments and a lemma on the eigenvalues. 
Lemma B.2 (Truncated moments [T71 Theorem VIII. 9. 2]). For every p > a, 

E[|Xi,iri{|x,,,Kt}] ~c(p)L(t)tP-" 
where c{p) ■.— a/{p — a). In particular, E i|^a„}] ^ c(p)afj/n. 

Lemma B.3 (Schatten bound [291 proof of Theorem 3.32]). If A is an n x n complex Hermitian 
matrix then for every < r ^ 2, 

n n In \ 

^|Afc(A)r<^K]|A,^,f • (B.l) 
k=l i=l \j=l J 

Proof of Lemma \B.1\ 
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Proof of (i). Let us fix r > 0. By definition of fix we have 

/•CO -1 ^ 

•^0 k=l 



By using (jB.ll) we get for any ^ r ^ 2, 

/ W ^^a-'xidt)^Z,,■.= -^Yn,^ where Y^,r.= \^a-^\X,^j\^ \ . 

We need to show that (2'„)„^i is a.s. bounded. Assume for the moment that 

supE(r„''i) < oo (B.2) 

for some choice of r. Since l^n,i, . . . , Yn^n are i.i.d. for every rt ^ 1, we get from (|B.2p that 
E((Z„ - EZ„)4) = n-^E 1^ y„,, - Ey„,,^ j = 0(n-2). 

Therefore, by the monotone convergence theorem, we get ^{J2n>i(^'n- ~ lEZn)*) < oo, which gives 
J2n^ii^n — ^Zn)'^ < OO a.s. and thus Z„ — EZ„ — > a.s. Now the sequence (EZ„)„^i = (EF„_i)„^i 
is bounded by (IB. 21) and it foUows that (Zn)n^i is a.s. bounded. 

It remains to show that (jB.2l) holds, say if < 4r < a. To this end, let us define 

n 

Sn,a,b := Xl"«^l^i'jl^^{a^'|Xi..P6[a,&)} ^^^^ry a < b. 
Now y;f_i (5'„,o.oo)^'' = {Sn,o,i + S'„4,oo)^'' and thus, 

nri^) ^ 2^^-'{nsi:o,i)+nsi:,,^)}- (b.s) 

We have sup„ E(5^'q i) < oo. Indeed, since 2r < 1, from the Jensen inequality, 

nsi:o,i) ^ (E5„^o,i)"^ 

and, by Lemma IR21 ES'„,o,i ~n a/(2 — a). 

To deal with the second term of the right hand side of (|B.3p . we define 

M„ := max a'^ |1 , -i,^ and iV„ #{1 sC j < n s.t. a-^\Xi,j\ > 1}. 

From the Holder inequality, if l/p-\- l/q= 1, we have 

IE(d,oo) ^ E (^fiVe) ^ (EiV^'-P)^/^ (EM,f «)'/^ . (B.4) 

Recall that P(|Xi_2| > = (1 + ^ 2/n for large enough n. Using the union bound, for 

large enough n, 

In particular for any 77 > 0, sup„ EA^^ < 00. Similarly, since L is slowly varying, for large enough 
n and alH ^ 1, 

P(M„ ^ t) < nPdAi^zl > tan) = na,7"i""L(a„t) < 2i"". 



It follows that if 7 < a, sup„ EM^ < 00. Taking p and (7 so that 4rg < a, we thus conclude from 
(lEl that sup„E(52r.^) <oo. 
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Proof of (ii). Recall that for any a e [1,2), k„ = w„na„^. Then, by using (IB. II) we get for any 
^ r ^ 2, 

r/2 

|2 



/ IC f^.„K{dt) <; - V r„,, where y„., := V a-'\X,. 

Jo n \ p, ) y^.^^ 

From (j2.18p (for 1 < a < 2) and p.6p (for a — 1), there exists c > such that a.s., 



nmsup max ) < c. 



Hence for all n large enough, 



and we conclude by using the same argument as in the proof of (i). □ 
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Figure 1 . Histograms of scaled ESDs illustrating the convergence stated by The- 
orems O and O for the following values of a: 0.25, 0.50, 0.75, 1.00, 1.25, 1.50, 
1.75, 2.00. Here n = 5000 and C is the law of V~^/" where V is a uniform 
random variable on (0,1). The first three plots are the histogram of the spec- 
trum of a single realization of K. The fourth plot corresponds to a = 1 and 
is a histogram of the spectrum of a single realization of log(n)7^. The four last 
plots are the histogram of the spectrum of a single realization of ■ In order 
to avoid scaling problems, an asymptotically negligible portion of the spectrum 
edge was discarded: only A[iog(„)j , • ■ ■ , A[„_iog(„)j were used. The liability of these 
simulations is questionable due to the lack of numerical precision for heavy tails. 



SPECTRUM OF MARKOV CHAINS ON COMPLETE GRAPH WITH HEAVY TAILS 



(Ch. Bordenave) IMT UMR5219 CNRS, Universite de Toulouse III, France 
E-mail address: Charles. bordenave(at)math. univ-toulouse.fr 
URL: http : //www . math . univ-toulouse . f r/"bordenave/ 

(P. Caputo) DiPARTIMENTO Di Matematica, Universita Roma Tre, Italy 
E-mail address: caputo(at)mat. uniroma3.it 
URL: http : //www .mat . uniromaS . it/users/ caputo/ 

(D. Chafai) LAMA UMR8050 CNRS, Universite Paris-Est Marne-la-Vallee, France 
E-mail address: dj alii (at) chafai. net 
URL: http://djalil.chafai.net/ 



